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One of the surprising results from the Hubble Space Telescope was the 
discovery that many of the most massive galaxies at redshift z = 2 are very 


compact, having a half-light radius of only 1-2 kpc. The interpretation is that 
massive galaxies formed inside out, with their cores largely in place by z = 2 
and approximately half of their present-day mass added later through 
minor mergers. Here we present a compact, massive, quiescent galaxy 

at a photometric redshift of Zphor = 1.947? 73 with a complete Einstein ring. 


—0.17 


The ring was found in the James Webb Space Telescope COSMOS-Web survey 
andis produced by a background galaxy at Zphot = 2.987042. Its 1.54” 


—0.47° 


diameter provides a direct measurement of the mass of the ‘pristine’ core 
ofa massive galaxy, observed before the mixing and dilution ofits stellar 


population during the 10 Gyr of galaxy evolution between z = 2and z = 0. We 
find a mass for the lens Miens = 6.5+?7 x 10" M, within a radius of 6.6 kpc. The 
stellar mass within the same radius is Mstars = 1.1703 x 10" M, for a Chabrier 
initial mass function and the fiducial dark matter massis Mam = 2a' x 10" M,. 
Additional mass appears to be needed to explain the lensing results, either 

in the form of ahigher-than-expected dark matter density or abottom-heavy 


initial mass function. 


The galaxy and its ring were identified in James Webb Space Telescope 
(JWST) NIRCam observations in the context of the COSMOS-Web 
project’, a public wide-area survey using the F115W, F1SOW, F277W 
and F444W filters. A visual inspection of a mosaic generated from 
the F115W, F277W and F444W data available as of 15 April 2023, cover- 
ing atotal area of 0.35 deg’, readily revealed the object (Methods). The 
NIRCam images containing the galaxy were resampled’ to a common 
0.025” per pixel grid for analysis. 

The object, dubbed JWST-ERI, is shown in Fig. 1a. It consists of a 
compact early-type galaxy (JWST-ER1g) and a complete Einstein ring 
(JWST-ERIr) with two conspicuous red concentrations. The lensed 


galaxy probably has a red centre and a blue disk, with parts of the disk 
producing the ring. The diameter of the centre of the ring is 1.54” + 0.02”. 
JWST-ER1 joins a large number of known Einstein rings*", although most 
arenot complete. Like other strong lensing configurations, Einstein rings 
canbeusedto reconstruct high-resolution images of lensed background 
galaxies, using ray-tracing techniques’. However, the unique value of 
Einstein rings is what they tell us about the lenses themselves: given 
the redshifts of the lens and source, they provide a model-independent 
measurement of the enclosed mass within the radius of the ring’. 

We obtained five-band photometry of the lens by fitting it with a 
Sérsic model’, masking the ring and keeping the structural parameters 
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Fig. 1| A complete Einstein ring identified inJWST images. a, Colour image of JWST-ERI, created from the NIRCam F115W, F150W and F277W data. b, Model of the 
galaxy, with an effective radius of r, = 1.9 kpc. c, Residual of the fit. Each panel spans 4.1” x 4.1”. The coordinates of the lens are RA = 10 h 0O min 24.11s, 


dec. = 1° 53’ 34.9’ (J2000). 


fixed in all bands. The effective radius of the galaxy r, = 0.22” + 0.02” 
and its Sérsic index n=5.0 + 0.6. The total magnitudes of the galaxy 
are given in Table 1 and the spectral energy distribution (SED) is 
shown in Fig. 2a. There is a pronounced break between the F814W and 
F115W bands, leading to a well-constrained photometric redshift of 


z= 1.9403 for the lens (Methods). The photometric redshift 


exceeds the spectroscopic redshift of the most distant known lens, a 
z= 1.62 galaxy inacluster®. The source redshift is less well constrained. 
We split the source into two photometric masks, one containing 
the blue ring and one covering both of the red knots. The blue ring 
shows no strong features and has a redshift of Zpnot = 2807) The 
SED of the red knots has a clear break between F150W and F277W, 


and a better-constrained redshift of Zpnot = 2.9819 4? (Fig. 2b). 


The lensing galaxy appears to be a textbook example of the class 
of massive quiescent galaxies at z = 2. Its rest-frame colours, U-V = 2.10 
and V-J = 1.3, place it comfortably in the quiescent region of the 
z= 2 UVJ diagram’. The best-fit stellar population parameters from 
the Prospector” fit imply an age of 1g Gyr and alow star formation 
rate of 4*? M, yr™. The Prospector total stellar mass of JWST-ER1g is 


1.3? 3 x 10" M, for a Chabrier" initial mass function (IMF), and its 
observed effective radius corresponds tor, = 1.9 + 0.2 kpc. This makes 
the galaxy quite compact, just like other quiescent galaxies at these 
redshifts», and it falls on the canonical size-mass relation of 
quiescent galaxies“. The galaxy is almost perfectly round and 
there areno obvious star-forming regions, tidal tails or other irregulari- 
ties in the residuals from the GALFIT fit. 

We now turn tothe mass of JWST-ER1g as inferred from the radius 
of the Einstein ring. The photometric redshifts of the lens and source, 
combined with the radius of the Einstein ring, give a total mass of 
Mens = 6.5732 x 10" M,within r= 6.6 kpc (Methods). The stellar mass 
within the Einstein radius is 0.79 times the total mass as determined 
by GALFIT and Prospector, that is, Mars = (11703 x 10") M, fora 
Chabrier IMF. There is a large difference between the lens mass and 
the Chabrier stellar mass of JWST-ER1g, with the lens mass a factor of 
sot higher than the stellar mass. This is the central result of our 
study (besides the report of the discovery of JWST-ER1), and in the fol- 
lowing we discuss several possible contributors to the lensing mass. 

It is unlikely that a significant fraction of the lensing mass is in 
the form of gas. Observations of lensed quiescent galaxies”, as well as 
simulations'®”, have consistently found low gas masses (<10 Mo) for 
massive quiescent galaxies at these redshifts. Furthermore, a total gas 
mass of 3 x 10" M, within 6.6 kpc corresponds to sucha high projected 
gas density that a high star formation rate is inevitable. The average pro- 
jected surface density would be -2,200 M, pc”, and according to the 
Kennicutt-Schmidt relation” the corresponding star-formation-rate 


Table 1 | Structural parameters of the lens 


Filter F, (pixels) n b/a PA 
F15W 7.9+0.7 41+0.3 0.94 77 
F150W 9.9+0.5 4.9+0.2 0.96 -15 
F277W 8.8+0.4 5.3+0.2 0.98 -16 
F444W 8.6+0.4 5.6+0.2 0.99 23} 


nis the Sersic index. PA is the position angle. 


surface density is X = 14 M, yr’ kpc”. The total star formation rate 
within the ring would be ~2,000 M, yr“, which is three orders of 
magnitude higher than derived from the Prospector fits and 30 times 
higher than an upper limit derived from Spitzer/MIPS 24 um data 
(Methods). This is a rough estimate, with the actual star formation 
rate depending on the distribution and temperature of the gas, but 
the point is that J)WST-ER1g would not be quiescent but rather a strong 
starburst galaxy. 

There is of course dark matter within the Einstein ring, and 
with standard assumptions this explains about half of the difference 
between the lensing mass and the stellar mass. Assuming a 
Navarro-Frenk-White (NFW) profile” and the stellar mass-halo mass 
relation” for z= 2, the dark matter mass within the Einstein radius 


is Mam = 2.675% x 10" M, (Methods). As shown in Fig. 3, this leaves 


2.8+34 x 10" M, unaccounted for. An explanation for this mild 


discrepancy is that the dark matter density within the Einstein radius 
is a factor of -2 higher than expected from scaling relations. The ‘extra’ 
dark matter can come in two forms. First, the total halo mass could be 
higher than what is indicated by the canonical stellar mass—halo mass 
relation. Asecond option is that baryonic processes have led to a dark 
matter profile that deviates from the NFW form. The final profile can 
be steeper or shallower in the central regions, depending onthe balance 
between cooling and feedback”. 

Looking closer, both options are somewhat unlikely in the specific 
case of JWST-ER1g. As detailed in the Methods section, the total halo 
mass would have to be very high, close to Maio 10“ Mo, and only a few 
halos of that mass are expected to exist in the surveyed volume. Turn- 
ing to baryonic processes, they tend to alter the dark matter profile on 
the spatial scales where the baryons are located: specifically, signifi- 
cantly steeper profiles are expected in regions where the stellar mass 
dominates”, thatis, at radii of at most r,. The dark matter mass within 
1.9 kpc is only 3 x 10° M, for a 10” M, halo with a NFW profile, and 
even if this were enhanced by a factor of 2-3, it would not be enough to 
account for the missing mass within the Einstein radius. 
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Fig. 2| Photometry of the lens and source. Top panels: images in the HST/ACS 
F814W band and the JWST/NIRCam F115W, F150W, F277W and F444W bands. They 
are shown at a common F; scale. a, SED of the lens galaxy, determined from forced 
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and a total stellar mass of 1. ie -3 x 10" M, (for a Chabrier IMF). b, SED of the 
lens galaxy, with the summed flux of the two red knots shown in black circles and 
the blue ring in grey circles. The red knots provide a reasonably well-constrained 
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Fig. 3 | Comparison of stellar mass to lensing mass for different IMFs. a, The 
three IMFs that are considered in this study: a Chabrier IMF", a Salpeter IMF” 
and an IMF that is steeper than Salpeter with a slope of -2.7. Broken lines are for 
stars above the turn-off mass. b, Comparison of the stellar mass to the lensing 


Lensing mass, inside r, (M,) 


mass (open symbols) and the lensing mass minus the fiducial amount of dark 
matter (solid symbols), for the three IMFs. The dashed line indicates a one-to-one 
relation. Data points are measurements + s.d. rẹ, Einstein radius; x,, number of 
stars; m, mass of stars; M,,,, total mass. 


Anintriguing alternative is that the missing mass is in the form of 
low mass stars, and that the stellar IMF needs to be adjusted: stars with 
masses M = 0.5 M, and below dominate the total mass but contribute 
less than 5% to the light”. Rather than simply scaling the mass, we 
refit the photometry in Prospector with two bottom-heavy IMFs: 
the Salpeter form”, with a slope of -2.3 and no turnover, and a 
‘super-Salpeter’ IMF with a slope of -2.7. These IMFs are illustrated in 
Fig. 3a. We note that these parameterizations are not unique, as the low 
mass slope is degenerate with the low mass cut-off. Furthermore, 
top-heavy IMFs can lead to high M/L ratios too if the mass is dominated 


by stellar remnants, although even for very flat IMFs this only occurs 
at ages greater than 3 x 10° Gyr (ref. 28). With these caveats in mind, we 
find that the stellar mass within the Einstein radius is 2. One > x 10"M, 
for a Salpeter IMF and 4.0+°¢ x 10" M, for the super- Salpeter IMF. 
As shown in Fig. 3b a model that combines a super-Salpeter IMF with 
a standard dark matter halo matches the lensing mass exactly, with a 
Salpeter IMF also providing a good fit. 

The probable descendants of compact quiescent galaxies at z = 2 
are massive early-type galaxies” >”, and the central regions of these 
galaxies may indeed have IMFs that are more bottom heavy than the 
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Chabrier IMF. The evidence largely comes from gravity-sensitive 
absorption lines”, kinematics™ and gravitational lensing”. Out 
side the central regions there appears to be a gradual transition toa 
Chabrier IMF*®”, as expected if a significant fraction of the mass in 
the outskirts was accreted through minor mergers. Quantitatively, the 
excess stellar mass compared to a Chabrier IMF reaches a factor of -3 
inthe centres of massive galaxies, with a power-law slope of -2.7 found 
for the galaxy NGC 1407 from a detailed non-parametric analysis*°. 
Super-Salpeter slopes of -2.7 have also been proposed on theoretical 
grounds“. We infer that a steep IMF for JWST-ER1 would be consistent 
with estimates in the central regions of early-type galaxies, particularly 
when mixing and dilution due to mergers and projection effects are 
taken into account. 

While this consistency is encouraging, IMF measurements are 
difficult and often indirect, and the question of IMF variation in the 
central regions of elliptical galaxies is still debated*’*’. Furthermore, 
and of direct relevance to JWST-ER1g, bottom-heavy IMFs are in 
some tension with comparisons of dynamical masses to stellar masses 
of z= 2 galaxies, which tend to prefer bottom-light IMFs such as 
the Chabrier form***’. On the other hand, our results are qualita- 
tively consistent with the most similar system to JWST-ER1, which is a 
z= 1.525 lens and partial Einstein ring that is best fit with a Salpeter IMF°. 

The combination of lensing with kinematics can break some of the 
degeneracies between the dark matter profile and the stellar mass, as 
has been demonstrated at lower redshifts“. This should work particu- 
larly well for JWST-ER1 as the effective radius of the galaxy is a factor 
of 3.5 smaller than the Einstein radius. Future NIRSpec observations of 
JWST-ER1 could provide the velocity dispersion of the galaxy, and pin 
down the redshifts of the lens and source. 


Methods 

Discovery 

JWST-ER1 is located in the COSMOS-Web JWST data’, as described 
in the main text. We reduced and aligned the NIRCam images witha 
software pipeline that was previously developed for the Hubble Space 
Telescope (HST) imaging and was modified for the JWST instruments”. 
Existing HST/ACS F814W imaging from the original COSMOS project*® 
and datasets at other wavelengths were processed in the same way, 
so that all space-based images are aligned to a common astrometric 
frame. The galaxy was found from a visual inspection of a mosaic that 
was generated from the F115W, F277W and F444W data. 

It is not the first Einstein ring that was found in the COSMOS 
field; there are at least two others, along with several more candidates”. 
This raises the question of why JWST-ER1 had not been noticed before. 
The main reason is that both the source and the lens are faint in the 
optical, and that existing HST data in the near-infrared—while show- 
ing the lens—are not deep enough to show the source. In Supplemen- 
tary Fig. 1 the pre-JWST high-resolution data are shown: the HST/ACS 
F814W fromthe original COSMOS programme‘ and a short-exposure 
HST/WFC3 F160W image from 3D-DASH, a wide-field survey with the 
drift-and-shift (DASH) technique’. With the benefit of hindsight, the 
characteristics of an Einstein ring can be glimpsed: a compact red 
galaxy near the centre of a blue ring. 


Is it alens? 
We consider the possibility that the system is not a gravitational lens but 
aring galaxy, suchas Hoag’s object”. Star-forming rings can be created 
in head-on collisions” and there is a small galaxy to the southwest of the 
ring that could be the perturber. The most obvious argument in favour 
of the lensing interpretation is that the photometric redshift of the ring 
is higher than that of the central galaxy (see main text). However, the 
redshift of the ring is uncertain, and it might be possible to fit both the 
lens and the ring with a model at z = 2.1. 

Here we highlight the morphology of the ring. In Supplementary 
Fig. 2 we show an enlarged, high-contrast colour image generated from 


the FISOW and F444W data, after subtracting the best-fitting model for 
the central galaxy. There are several symmetries in the image: as well 
as the two bright red knots it appears that two blue knots are also multi- 
ple imaged. The most compelling argument for lensing is the morpho- 
logy of the red knots (presumably the bulge of the lensed galaxy): they 
are stretched into mirrored arcs on each side of the galaxy, something 
that cannot be explained in collisional ring scenarios. 


Structural parameters 

We fit the lens galaxy with the GALFIT code” to determine its structure 
and in preparation for measuring its photometry. We use cutouts of 
4.1” 4.1” with 0.025” per pixel sampling in the NIRCam bands and 
0.05” per pixel sampling in the ACS F814W band. The presence of the 
ring makes it difficult to measure the size, Sérsic index and background 
level simultaneously. We therefore first measure the background level 
in each band from the outer edge of the cutout, iteratively rejecting 
outlying pixels, and subtract this value. Next a mask is created for the 
ring, by selecting pixels in the ring area above a flux threshold and then 
expanding the mask using a5 x 5 pixel kernel. 

The fit is performed on the F115W, FISOW, F277W and F444W 
images (the signal-to-noise ratio in the F814W image is too low 
for a stable fit). Free parameters are the position, Sérsic index, effec- 
tive radius, total magnitude, axis ratio and position angle. We use 
the WebbPSF tool (https://www.stsci.edu/jwst/science-planning/ 
proposal-planning-toolbox/psf-simulation-tool) to create point 
spread functions (PSFs) for each filter and position. We verified that a 
well-exposed nearby star does not lead to qualitatively different results. 

The structural parameters are listed in Table 1. The parameters 
inthe four bands are in good agreement, despite the factor of 4 range 
in wavelength and resolution going from F115W to F444W. The aver- 
age effective radius r, = 8.8 + 0.8 pixels, or 0.22” + 0.02”, where the 
root mean square of the four individual measurements is taken as the 
uncertainty. The Sérsic index n = 5.0 + 0.6. The axis ratio is very close to 
land there is no consistent position angle between the bands; in what 
follows we therefore assume that the axis ratio b/a =1.0. 


Photometry 

Total magnitudes of the lens are determined by fitting the five bands 
(now including ACS F814W) with GALFIT, holding all parameters except 
the total magnitude fixed to the average values determined above. This 
constrained (or forced) fit ensures that the relative fluxes between the 
bands are measured in a self-consistent way, and not compromised by 
PSF or aperture effects. The results are listed in Supplementary Table 
1, with 0.05 mag systematic error added in quadrature to the random 
errors. For the comparison of the lensing mass to the stellar mass it is 
not the total flux but the projected flux within the Einstein radius that 
matters. Using a model profile that is not convolved with the PSF we 
determine that 79% of the total flux is within the Einstein radius. For 
convenience the magnitudes within the Einstein radius are listed ina 
separate column. We tested that simple aperture photometry on the 
galaxy, with the ring masked, gives a redshift and M/L ratio that are 
within the uncertainties of the fiducial values. 

Photometry of the ring is performed by simply summing the flux 
in apertures. Two apertures are used: one covering both of the red 
concentrations within the ring and one covering the rest of the ring. 
Noattemptis made to correct for the PSF variation between bands, but 
the apertures are purposefully made large enough to mitigate these 
effects. We use the photometry of the ring to derive an approximate 
redshift, and we caution against using it to determine detailed stellar 
population parameters of the lensed galaxy. The magnitudes for the 
two apertures are listed in Supplementary Table 2. 


Prospector fits 
The redshift of the lens and its stellar population parameters are deter- 
mined jointly using the Prospector inference framework”, specifically 
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the Prospector-a@ model™ and the MIST stellar isochrones**”° from 


Flexible Stellar Population Synthesis”. Prospector-a@ describes the 
star formation history (SFH) non-parametrically via mass formed in 
seven logarithmically spaced time bins, and assumes a continuity 
prior to ensure smooth transitions between bins**. We additionally 
adopt a dynamic SFH(M, z) prior” that follows the observed cosmic 
star-formation-rate density, favouring rising SFHs in the early universe 
and falling SFHs in the late universe, with a mass-based adjustment to 
reflect downsizing. The model consists of 18 free parameters, including 
the form of the attenuation curve, and sampling is performed using 
the dynamic nested sampler dynesty°’. The parameters for the lens 
are determined from the photometry inside the ring. We report the 
posterior median of the inferred physical parameters in Supplementary 
Table 3, assuming a Chabrier IMF. The uncertainties reflect the 16th 
and 84th percentiles. 

The uncertainties in the redshift and mass may seem suspiciously 
small given that we only have five photometric datapoints. The reason 
why the key parameters are so well constrained is that the photometry 
tells us only one thing, but it does so precisely: there is a large breakin 
the SED at 1.2 um. The constraints on the redshift and M/L ratio follow 
directly from this. We performed two robustness tests to determine 
howsensitive the results are to the specifics of our methodology. First, 
removing the SFH prior leads to negligible differences to the redshift 
and mass, and all the posterior medians are consistent within lo. The 
only notable change is that the prior decreases the uncertainty on the 
star formation rate. This behaviour is expected: at this redshift 
and mass the prior prefers a falling SFH, consistent with the observed 
high mass (that is, high previous star formation rate) and low 
current star formation rate. Second, determining the redshift with 
the EAZY code“ (which uses a pre-rendered set of templates) gives 
Zphot = 1.914?-18 and no viable secondary solutions, in good agreement 
with our fiducial value. 

The lensed galaxy is modelled in the same way as the lens, except 
that the scale-dependent SFH prior is not included owing to the lensing 
magnification. The main goal is to determine the redshift of the lensed 
galaxy. For completeness we list stellar population parameters for the 
two apertures on the ring as well in Supplementary Table 3, although 
they are not used in the analysis. 


Obscured star formation 

The low star formation rate of JWST-ER1g derived above implies alow 
gas surface density, and hence a low contribution of gas to the total 
mass budget within the Einstein ring. However, the Prospector fits do 
not provide strong constraints on the amount of star formation that is 
optically thick. The field has been observed with Spitzer/MIPS, as part 
of theS-COSMOS survey”, and we use the 24 um data to assess whether 
JWST-ER1g has a hidden obscured star burst. 

The S-COSMOS 24 um image is shown in Supplementary Fig. 3. The 
galaxy is not detected. We determine an upper limit to the star forma- 
tion rate from a redshift-dependent relation between observed 24 ym 
flux and total infrared luminosity that was calibrated with Herschel 
data®°*. The 3o upper limit is 63 M, yr”. 


Comparison to other z = 2 galaxies 

As noted in the main text, JWST-ER1g is a typical example of the class 
of massive, quiescent z= 2 galaxies. This is demonstrated explicitly 
in Supplementary Fig. 4. Supplementary Fig. 4a shows that the gal- 
axy falls in the quiescent region of the UV) diagram. The boundaries 
are the averages of the z=1.75 and z= 2.25 limits determined for the 
NEWFIRM Medium Band Survey’. Itis relatively red within the quiescent 
region, indicating an old age and/or some dust, as also implied by the 
Prospector fit. In Supplementary Fig. 4b the galaxy’s size is compared 
to the canonical size-mass relations” for quiescent and star-forming 
galaxies, again taking the average of the listed relations for z= 1.75 and 
z= 2.25. The galaxy falls on the relation for quiescent galaxies. 


Lensing mass 
The mass within the Einstein radius is given by: 


6??D,D, 


MG leary a 
iS 


(1) 


with 6 the observed Einstein radius in radians, c the speed of light, 
D, the angular diameter distance to the lens, D, the angular diameter 
distance to the source and Gthe gravitational constant. The parameter 
D,, is the distance between the lens and the source, which is: 

14+z, 


LD, 2) 


D= Peg 
S 


with z the redshift of the lens and z, the redshift of the source in a flat 
Universe“. The uncertainties are determined numerically, by drawing 
values of z,,z,and @ from their probability distributions and calculating 
M(<6) for each set of draws. 

The high lens mass is driven by the large diameter of the Einstein 
ring combined with the relatively high redshift of the lens. Forcing 
Z,=1.5 (which is outside of the full posterior distribution of 5,000 
samples) lowers the mass to M,,,; = 4.1 x 10" M,, but also lowers the 
derived Chabrier stellar mass to M,,,,,= 0.6 x 10" M,. The ratio of the 
lensing mass to the Chabrier mass is ~7, which is very similar to the 
results for z= 1.94. 

The source redshift is the most uncertain parameter in equa- 
tion (1). The lensing mass is lower for higher source redshifts, but is 
3.7 x 10" M, even for z, = 5. The uncertainty in the source redshift also 
causes an asymmetry inthe error distribution of Mens with a tail to very 
high masses. This is because the mass increases rapidly when z, = z: the 
mass is >10” M, if z, < 2.5, and reaches 4 x 10” M, for z, = 2.1. 


Dark matter contribution 

The projected dark matter mass within the ring can be calculated 
by integrating a NFW profile” along acylinder witha radius of 6.6 kpc 
(ref. 66). The scaling log c = 0.81 — 0.09(log M,i, — 12s used to determine 
the concentration as a function of halo mass, with M,,, the virial 
mass”. The resulting relation between projected dark matter mass 
within the ring and total halo mass is shown in Supplementary Fig. 5. 

The relation is shallow, owing to the decreasing concentration 
with halo mass. We estimate the dark matter contribution to the 
lensing mass from the halo mass-stellar mass relation”. We find 
Movo = 1.072 x 108 Mo, where Moo is the mass within the radius where 
the overdensity is a factor of >200, with the relatively large uncertainty 
driven by the steepness of the relation inthis regime. The correspond- 
ing projected dark matter mass within 6.6 kpc is Mam = 2.6*G° x 10"M, 
for a NFW halo. 

The solid horizontal line indicates the difference between the 
lensing mass and the stellar mass of JWST-ER1g for a Chabrier IMF. To 
explain the missing mass entirely with dark matter, the NFW halo mass 
would have to be -7 x 10” M,. Halos of this mass at z = 2 are progenitors 
of clusters at z= 0. The number density of halos with Mj, > 7 x 10° M, 
at z=1.94 is 2 x 107” h? Mpc”, corresponding to 1.4 in the redshift 
range 1.75 < z < 2.25 in the 0.35 deg’ of the available COSMOS-Web 
area’. Halos with slightly lower masses are of course more common, 
and still consistent with the lensing constraints. The lower 1g bound 
on the lensing mass corresponds to a halo mass of Moo > 3 x 10° M, 
(Supplementary Fig. 5), and there are ~15 such halos in the COSMOS- 
Web area. 


Environment of JWST-ER1 

Gravitational lensing is sensitive to the weighted integral of all mass 
between the source and the observer, and we briefly consider whether 
nearby galaxies or structures along the line of sight could contribute 
to the mass. We also consider whether JWST-ER1g is the central galaxy 
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of the progenitor ofa cluster (see above). The immediate environment 
of JWST-ER1 is shown in Supplementary Fig. 6, as generated from the 
NIRCam F115W, F277W and F444W images. The region does not stand 
out in any way; the galaxy is either isolated or in a sparse group, but 
notin a massive cluster. Furthermore, there are no other bright galaxies 
projected along the line of sight. We infer that the contributions from 
other galaxies to the 6.7 x 10" M, mass within the Einstein radius are 
almost certainly negligible. 


Data availability 
The COSMOS-Web data are publicly available from the STScl MAST 
Archive. 


Code availability 
We have made use of standard data analysis software in the Python and 
IRAF environments, and the publicly available code GALFIT™. 
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